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^ Abstract 



We compare the eigenvalues of the Dirac and Laplace operator on a two- 
dimensional torus with respect to the trivial spin structure. In particular, we 
' compute their variation up to order 4 upon deformation of the flat metric, study 

fH , the corresponding Hamiltonian and discuss several families of examples. 
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1 Introduction 



We consider a two-dimensional torus T 2 equipped with a flat metric g a as well as a 
conformally equivalent metric g, 



9 = h 4 g Q . 

X 

Denote by A g the Laplace operator acting on functions and let D g be the Dirac 
operator. The following estimates for the first positive eigenvalues fJ-i(g) and A 2 (g) 
of A g and -D 2 are known (see [2] and [6]): 

\ A 2 (g ) 2 , A 2 (ff D ) /XjXgg) , . Ml(ffo) 

a T2 — < Ms) < -ra — . ~a — ^ miw < 



""max ""mm 'max ""mm 

where /i m i n (/i ma x) denotes the minimum (maximum) of the conformal factor. 



167T 

b.) m(g) < 



vol(T 2 , ff )' 
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In case the spin structure of the torus T 2 is nontrivial, the Dirac operator has no 
kernel and, moreover, there exists a constant C depending on the conformal structure 
fixed on T 2 and on the spin structure such that 

(see [7]). However, explicit formulas for the constants are not known. In this respect, 
the situation on T 2 clearly differs from the case of the two-dimensional sphere S 2 , 
where 

k 2, * - 47r 



K{g) > 



vol(S 2 ,<?) 
holds for any metric g (see [3], [7]). 

In this paper we compare Hi(g) and X 2 (g) for the trivial spin structure and a metric 
with S^-symmetry. We will construct deformations gE of the flat metric such that 
vo1(<7e) = vol (g ) and H\{gE) < A 2 (<?e) holds for any parameter E ^ near zero. 
For this purpose we calculate, in complete generality, the formulas for the first and 
second variation of the spectral functions fj,\ and A 2 . It turns out that for any local 
deformation of the flat metric, the first minimal eigenvalue of the Laplace opera- 
tor is always smaller than the corresponding eigenvalue of the Dirac operator up to 
second order. The question whether or not there exists a Riemannian metric g on 
the two-dimensional torus such that A 2 (g) < fJ^i(g) remains open. Denote by A 2 (g; I) 
and Hi(g; k) the first eigenvalue of the Dirac and Laplace operator, respectively, such 
that its eigenspace contains the representation of weight I respectively k. The dis- 
cussion in the final part of this paper suggests the conjecture that for same index 
/ = k, the eigenvalues X 2 (g,l) and fj,i(g,l), are closely related, more precisely, that 
the Laplace eigenvalue is always smaller than the Dirac eigenvalue and that their 
difference should be measurable by some other geometric quantity. 

The Dirac equation for eigenspinors of index I = can be integrated explicitly. In 
case of index / ^ 0, the Hamiltonian describing the Dirac equation is a positive 
Sturm-Liouville operator. First, this observation yields an upper bound for \ 2 {g,l). 
On the other hand, it proves the existence of many eigenspinors without zeros. 

We furthermore apply the general variation formulas in order to study the eigenvalues 
of the Laplace and Dirac operator for the family of metrics 

g E = (1 + Ecos(2irNt)) (dt 2 + dy 2 ) 



in more detail. In case N = 2, the Laplace equation is reduced to the classical 
Mathieu equation. A similar reduction of the Dirac equation yields a special Sturm- 
Liouville equation whose solutions we shall therefore call Mathieu spinors. We in- 
vestigate the eigenvalues of this equation and compute (for topological index / = 1) 
the first terms in the Fourier expansion of these Mathieu spinors. These computer 
calculations have been done by Heike Pahlisch and our grateful thanks are due to 
her for this. Furthermore, we thank M. Shubin for interesting discussions on Sturm- 
Liouville equations. 
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2 The first positive eigenvalue of the Dirac operator and 
Laplace operator 

Let g and g Q be two conformally equivalent metrics on T 2 . Then the Laplace and 
Dirac operators are related by the well-known formulas 

1 



A, 



h 4 



A , 



J_ grad(fr) 

where grad (h) denotes the gradient of the function h with respect to the metric g Q . 
Let us fix the trivial spin structure on T 2 . In this case the kernel of the operator D Q 
coincides with the space of all parallel spinor fields, in particular, any solution of the 
equation D (ip ) = has constant length. The kernel of the operator D g is given by 



ker(£> 9 



^Vo : A>(V>o) = 



The square D 2 of the Dirac operator preserves the decomposition S 



s+ e s- of 



the spinor bundle S. Moreover, D 2 acts on the space of all sections of with the 
same eigenvalues. Therefore, the first positive eigenvalue A 2 (g) of the operator D 2 
can be computed using sections in the bundle S + only. Any section G r(S' + ) is 
given by a function / and a parallel spinor field ip Q 6 

i> = (/ • m . 

The spinor field ip is £ 2 -orthogonal to the kernel of the operator D 2 if and only if 



J (V, \^o)dT 2 = |Vo| 2 J fh 4 dT 2 = 

J>2 ip2 



holds, where dT 2 and dT 2 = h A dT 2 are the volume forms of the metrics g D and g. 
The Rayleigh quotient for the operator D 2 is given by 

J \D 9 (i;)\ 2 dT 2 J |/ l -grad(/) + 2/grad(/ l )| 2 dr o 2 



T 2 



T 2 



T 2 



2 dT 2 



J fh 6 dT 2 

T 2 



Finally, we obtain the following formulas for the first positive eigenvalue Hi(g), A 2 (g) 
of the Laplace operator A 
the trivial spin structure: 



of the Laplace operator A g and the square D 2 of the Dirac operator with respect to 



|grad(/)| 2 dT 2 



Mi(ff) = inf < 



T 2 



21.4.^2 T 2 



J fh 4 dT 2 

T 2 



: / fh 4 dT 2 = 
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A?(<7) = inf 



|/i-grad(/)+2/-grad(/i)| 2 cff 2 



T 2 



T 2 



fh A dT 2 = 



2 h 6 dT 2 



T 2 



A direct calculation yields the formula 

|/ l -grad(/)+2/grad(/ l )| 2 dT 2 = J h 2 fA (f)dT 2 + J (4|grad (h)\ 2 '+\A (h 2 )) f 2 dT 2 . 



rp2 



T 2 



T 2 



Let us use this formula in case that / is an eigenfunction of the Laplace operator, 
i.e., 

A (f)=» 1 ( g )h 4 f. 

Then it implies the following inequality between the first eigenvalues of the Laplace 
and Dirac operator: 



A 2 (g) < fii(g) + 



J (4|grad(/ l )| 2 + lA (/ l 2 )) f 2 dT 2 

T 2 



If 

T 2 



2 h 6 dT 2 



We are now looking for L 2 -estimates in case the metric g admits an S 1 -symmetry. 
Indeed, let us suppose that the metric g is defined on [0, 1] x [0, 1] by 

g = h\t)g = h 4 (t)(dt 2 + dy 2 ), 

where the conformal factor h A depends on the variable t only. Moreover, we assume 
that the function h (t) has the symmetry 

h{t) = h(l-t). 

Then any function f(t) with f(t) = — /(l — t) satisfies the condition 

J fh 4 dT 2 = 

T 2 

and, consequently, yields upper bounds for ix\{g) and Ai(g): 



\f'(t)\ 2 dt 



J f 2 {t)h\t)dt 



B u L (g;f) 
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I (h(t)f'(t) + 2f{t)h'(t)fdt 

A?fo) < ° - == B u D (g;f) ■ 

J f\t)h\t)dt 



3 The first and second variation of fJ^i(g) and \\(g) 

We consider a Riemannian metric 

g = h\t) 9o = h\t){dt 2 + dy 2 ) 

on T 2 (0 < t < 1, < y < 1) and denote by E(pL\(g)) and E(X 2 (g)) the eigenspaces 
of the Laplace operator and the Dirac operator corresponding to the first positive 
eigenvalue. The isometry group S 1 acts on these eigenspaces and therefore they 
decompose into irreducible representations 

%i(s)) = E(ti)®-ffiE(y and £(A?( 5 )) = £(Zi)e---0£(i n ), 
where XX^) denotes the 1-dimensional 5 1 -representation of weight A;. 

Proposition 1: T/ie weights k 2 a of the first positive eigenvalue fii(g) of the Laplace 
operator are always bounded by one: 

k 2 a <l. 

The weights l 2 ^ of the first positive eigenvalue \ 2 (g) of the Dirac operator are bounded 
by one under the condition 



max 



h'(t) 



h(t) 



< t < 1 < 3tt. 



Proof: Suppose that 

f(t,y)=A(t)e 2 ^y 

is an eigenfunction of the Laplace operator, A g f = /j,\(g)f. Then the function A(t) 
is a solution of the Sturm-Liouville equation 

-A"(t) = {^(g)h\t)-A7r 2 k 2 a }A(t). 

Then consider the function 

F(t,y) = A(t)e 2m y 

and remark that 
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A 9 F = ^(g)F + 47r 2 (l-k 2 a )^. 



Since J FdT 2 = 0, we obtain, in case of the first positive eigenvalue, that 

->2 / I 7?|2^t->2 



A g (F)FdTg J \F\ 2 dT 2 

Mi < = Mi + 4^ 2 (1 - kl) — . 

J \F\ 2 dT 2 J \F\ 2 h 4 dT 2 

The latter inequality yields k 2 a < 1 immediately. The corresponding result for the 
Dirac operator follows from the formula 



A? < A? + — ((1 - ! 2 ) / \F\V; + '-A J |F(()| 2 ^ 



/ 



F\ 2 h A dT 2 K T2 



o 

T 2 



where we have already used the differential equation (**) for A that will be derived 
in the next paragraph. □ 

Solutions of the Laplace equation A g f = Hi(g)f are given by solutions of the Sturm- 
Liouville equation 

-A"(t) = {^(g)h 4 (t) - 4n 2 k 2 }A(t) (*) 

with k = 0, ±1. In a similar way we can reduce the Dirac equation to an ordinary 
differential equation. The Dirac operator D g acts on spinor fields via the formula 



h 2 {t) { i o ) * ^ h 3 {t) { i o y h 2 (t) { i o 



Suppose that a spinor field tp € r(S' + ) is a solution of the equation D 2 (tp) = X 2 (g)ip. 
Then tp is given by a solution of the Sturm-Liouville equation 

-a-v = (xjto^w + m^pg _ ^ + 4rf m j . 4(t) M 

with Z = 0, ±1. In case Z = 0, this equation can be solved. 

Proposition 2: XTie eigenvalues of the Sturm-Liouville equation (**) /or Z = are 
given by (n G Z) 
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x2 4vr 2 n 2 
A = 



2 • 

h\t)dt 



\0 



Proof: The Sturm-Liouville operator 

1 d 2 h{t)h"{t) - 2(h'(t)) 2 

~ 'hHt)^ hP{t) 

admits a square root, namely 

^<->-j4iw'>->- 



Since we have ^h 2 (t)A(t)A(t) = 0, any solution of the equation 

' d (h(t)A(t)) = XA(t) 



h 3 (t) dt 
satisfies the condition 

. . const 

^ = w 

Consequently, it makes sense to define a function / : K 1 — > M 1 by the formula 

h(t)A{t) = e if(t \ 

for which we easily obtain the differential equation f'(t) = Xh 2 (t). But since A(t) 
is a periodic solution, we have the condition 

l i 



2im = J f'(t)dt = \J h 2 (t)dt 



o o 
for some integer n € Z, thus yielding the result. 

Corollary: Let \ 2 {g) be an eigenvalue of the Dirac operator on the two-dimensional 
torus T 2 with respect to the trivial spin structure and a Riemannian metric 

g = h\t)(dt 2 + dy 2 ) 

with isometry group S 1 . Moreover, suppose that the eigenspinor is S 1 -invariant 
(1 = 0). Then 

X 2 (g)vo\ (T 2 ,g) > in 2 

holds. 
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1 

Proof: Since the volume is given by vol(T 2 ,<7) = J h 4 (t)dt, the inequality follows 



o 

1 



directly from the Cauchy-Schwarz inequality if h 2 (t)dt J < J h 4 (t)dt and the 

Vo Jo 

previous Proposition. □ 

Remark: This corollary should be compared with the following fact. Fix a re- 
presentation S(/c) and denote by Hi{g\ k) the first eigenvalue of the Laplace operator 
such that its eigenspace contains the representation £(&). In case k ^ the solution 
A(t) of equation (*) is positive (see [8], page 207) and consequently the inequality 



J (4irk 2 - m(g; k)h A {t))dt > 





is valid. We thus obtain the estimate 

4vr 2 fc 2 > / u 1 ( 5 ;fc)vol (T 2 ,g) 
and equality holds if and only if the metric is flat. In particular (k = ±1) we have 

4vr 2 > ^1 (5) vol {T 2 ,g) 

for the first positive eigenvalue of the Laplace operator in case that its eigenspace 
contains the representation S(±l). 



Let us introduce the Hamiltonian operator Hi defined by the Sturm-Liouville equa- 
tion (**) for A 2 = 0: 



dt 2 h(t) h 2 (t) 

Proposition 3: For I 7^ 0, the Hamiltonian operators Hi are strictly positive. H Q is 
a non-negative operator. 

Proof: A direct calculation yields the formula 

2 



\h(t)J h(t) J V h(t) h(t) J 




where (p(t) is any periodic function. The equation 2irl(p(t) — <p'(t) = does not admit 
a periodic, non-trivial solution in case I / 0. Consequently, Hi is a strictly positive 
operator for I 7^ 0. 

Corollary: Fix an S 1 -representation Let -^i(s>0 ^ e ^ e fi rs ^ eigenvalue of 

the Dirac operator on the two-dimensional torus T 2 with respect to the trivial spin 
structure and an S 1 -invariant metric 
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g = h\t)(dt 2 + dy 2 ) 



such that the representation E(0 occurs in the decomposition of the eigenspace. Then 
the multiplicity o/E(0 * s one an d the eigenspinor does not vanish anywhere (I / 0). 

Proof: Since Hi is strictly positive, the eigenvalue A 2 (g) is the unique positive 
number A 2 such that 

inf spec(Hi - A 2 /i 4 ) = 0. 

The corresponding real solution of this Sturm-Liouville equation is unique and posi- 
tive (see [8], page 207). 

Corollary: For a fixed S l -representation E(0 denote by X 2 (g, I) the first eigenvalue 
of the Dirac operator such that the eigenspace E{\\{g, I)) contains the representation 
E(0- Then the inequality 



I 



1 (2nl<p(t)-p>(t)l dt 



h\t) 

\l(g,i)<° ■ 



J h 2 {t)^ 2 {t)dt 



holds for any periodic function f(t). 



Proof: Since inf spec(-£^ — X 2 (g, l)h A ) = 0, we have 
l l 



Hi (£) ^(t)dt - X 2 (g,l) J h 2 {t)^ 2 (t)dt > 



o o 
for any periodic function <p(t). 

In case of the flat metric g = dt 2 + dy 2 we have ^i(g ) = X\(g ) = Att 2 and 
E(»i(9o)) = E(0) © E(0) © E(l) © E(-l) = E{\ 2 { 9o )). 

The spaces correspond to the case that k = I = ±1 and are generated by 

the constant function. The two spaces E(0) are generated by the functions sin(27rf), 
cos(27rt). 

Notation: We introduce now a few notations which will be used throughout this 
article. Let us consider a deformation gE = h E (t)g Q of the flat metric g Q depending 
on some parameter E. We assume that 

h%(t) = h%(l-t) 
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holds for all parameters of the deformation. The eigenvalues Hi{g ) an d A 2 (<? ) of 
multiplicity four split into three eigenvalues 

fii(go)"{»i(E),» 2 (E),» 3 (E)} , X 2 (g )^{X 2 (E),X 2 (E),X 2 (E)}. 

The eigenvalue fJs(E) corresponds to the case that k = ±1, has multiplicity two, 
and its eigenfunction is a deformation of the constant function. The eigenvalues 
^i(E) / ^2(E) correspond to solutions of the Sturm-Liouville equation (*) and 
their eigenfunctions are deformations of sin(27rt) and cos(2-/rt), respectively. The 
situation is different for the Dirac equation: there, according to Proposition 2, the 
trivial S^-representation (I = 0) yields one eigenvalue X 2 (E) of multiplicity two and 
the non-trivial representations (7 = ±1) define in general two distinct eigenvalues 
\%{E), X\(E) of multiplicity one. However, in case = — t), the spectral 

functions X\(E) and X 2 (E) coincide. Obviously, for small values i?«0we have 

fi 1 (g E )=mm{ f i 1 (E),fi 2 (E),^(E)} , X\{g E ) = min {\\{E), \\{E), \\{E)}. 

We will compute the first and second variation of [i a (E) and A 2 (.E) at E = 0. For 
this purpose we introduce the following notation: Let A be a function depending both 
on E and t. Then A denotes the derivative with respect to E and A' the derivative 
with respect to t. Moreover, we expand the function h E (t) in the form 

h%(t) = 1 + EH(t) + E 2 G{t) + 0(E 3 ). 



Theorem 1: Consider a deformation 

g E = (l + EH(t) + E 2 G{t) + 0(E 3 ))g = h%{t)g 

of the flat metric on the torus T 2 such that h E (t) = h E (l — t). Moreover, suppose 
that for E / and k = the eigenvalues Hi(E) / fJ-2(E) are simple eigenvalues of 
the Sturm-Liouville equation (*). Then 

l l 

a. ) /^(O) = -8vr 2 J H{t) sin 2 (2nt)dt , /i 2 (0) = -8vr 2 J H{t) cos 2 (2nt)dt 

o o 
l 

/t 3 (0) = -4vr 2 y H{t)dt. 
o 

l 

b. ) A 2 (0) = A|(0) = A|(0) = -Ait 2 J H(t)dt. 

o 

In particular, we obtain 

Ai(0)+/i 2 (0) = 2/i 3 (0)=2A 2 (0). 
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Corollary: Suppose that the deformation 

g E = (l + EH(t) + E 2 G(t) + 0(E 3 ))g 
of the flat metric g Q satisfies the condition 

H{t) = H(l-t) 

as well as 



J H(t)sm 2 (2irt)dt^ J H(t) cos 2 (2irt)dt. 





Then, for all parameters E ^ near zero we have the strict inequality 

Mi (5b) < ^1(9 e) ■ 

Next we compute the second variation of our spectral functions under the assumption 
that the first variation is trivial. 

Theorem 2: Consider a deformation 

g E = (l + EH{t) + E 2 G{t) + 0(E 3 ))g = h%{t) 9o 
of the flat metric g Q on the torus T 2 and suppose that the conditions 

h%(t) = h%{l-t) 

and 

1 1 
J H{t) sin 2 (2irt)dt = J H(t) cos 2 (2irt)dt = 


are satisfied. Moreover, suppose that for E / and k = the eigenvalues H\(E) 7^ 
(J>2(E) are simple eigenvalues of the Sturm- Liouville equation (*). Then 

1 1 

a. ) /2i(0) = -16vr 2 J G(t) sin 2 (2vrt)^ - 16vr 2 J H{t)C{t) sin (2irt)db, 



where C{t) is the periodic solution of the differential equation 

C"(t) = -47r 2 H(t) sin (2vrt) - ^ 2 C{t). 

1 1 

b. ) fl 2 (0) = -167r 2 y G(t)cos 2 (2vrt) - 16vr 2 J H(t)C(t) cos(2Trt)dt, 



where C{t) is the periodic solution of the differential equation 
C"(t) = -4ir 2 H(t)cos(2irt) - 4vr 2 C(t). 
11 



c. ) /i 3 (0) = ~8tt 2 j G(t)dt - 8vr 2 J H{t)C{t)dt, 



where C{t) is the periodic solution of the differential equation 

C"(t) = -4ir 2 H(t). 

1 l 

d. ) A?(0) = -8n 2 J G{t)dt + 2ir 2 J H 2 (t)dt 



o o 



e.) A|(0) = A|(0) = -8vr 2 J G{t)dt + 4vr 2 ^ H 2 (t)dt - 8vr 2 ^ H{t)C(t)dt 



1 

-2vr y H'(t)C{t)dt, 
o 

where C{t) is the periodic solution of the differential equation 

C"(t) = -4ir 2 H(t)-irH'(t). 

Proof of Theorem 1 and Theorem 2: The formulas for the derivatives of \ 2 (E) 
are a direct consequence of Proposition 2. We will prove the variation formulas for 
A| and just remark that one can investigate the other spectral functions in a similar 
way. Moreover, since all the calculations we make are up to order two with respect 
to E, we may assume for simplicity that 

h%{t) = l + EH{t)+E 2 G{t). 

We compute 

h E h% - 2(h' E ) 2 _ 1_ E H" + EG" _ (H' + EG') 2 



h\ 4 (1 + EH + E 2 G) 16 (1 + EH + E 2 G) 2 

and, consequently, we obtain the formulas 

d ( h E h" E -2(h' E ) 2 \ 1 _d_ fh^\ = 1 

dE\ h\ ) E=Q 4 ' dE\h E J E=0 4 

d 2 (h E h' E -2(h' E ) 2 \ _ \ Q n _ H n H) _^ Hl) 2 



dE2 \ h l J E=0 2 

The spectral function X 2 (E) is defined by a periodic solution A E {t) of the Sturm- 
Liouville equation 

A%(Jt) = -\ 2 3 (E)h E (t)A E (t)- kEit ^ 
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with the initial conditions A§(0) = 4ir 2 ,A (t) = 1. Therefore, we obtain 

Al{t) = -A§(0) - 4vr 2 if(t) - 4vr 2 i (t) - ^H"{t) + 47r 2 i D (t) - Trff'ft) 
in this case and, consequently, 

1 

A 2 (0) = -4vr 2 y H(t)dt. 



Let us now compute the second variation in case that A 2 (0) = = A 2 (0). We 
differentiate the Sturm-Liouville equation twice at E = 0: 

%{t) = -A 2 (0) - 8vr 2 G(t) - ^ 2 H{t)A (t) + \{H'(t)) 2 - \(G"(t) - H"(t)H(t)) 
- l -H"{t)A {t)-2^H'(t)A (t)-A^ t (jL {ln{hEmE= ^ . 

Then we obtain 



ii i 
A 2 (0) = -8vr 2 J G(t)dt - 8tt 2 J H(t)A (t)dt + Q - ^ J (H'(t)) 2 dt 



1 1 
-- J H"(t)A (t)dt - 2ir J H'(t)A (t)dt. 



Since A (t) is a solution of the differential equation 



A>>{t) = -^ 2 H{t)-\H"{t)-itH>{t), 



we have 



i i i 

l -jH"{t)A {t)dt = -\J H(t)A'>(t)dt = +\j ' H{t)^H{t)+ l -H"{t)yt 



1 1 
= 2vr 2 J H 2 {t)dt - 1 J (H'(t)) 2 dt 



o 

and, consequently, we obtain 



ill l 
A 2 (0) = -8vr 2 J G(t)dt + 2tt 2 J H 2 (t)dt - 8vr 2 J H(t)A (t)dt - 2vr J H'(t)A (t)dt 



111 1 
= -8vr 2 J G{t)dt + 4vr 2 J H 2 {t)dt - 8vr 2 J H(t)C{t)dt - 2vr J H'(t)C(t)dt, 



13 



where C(t) := A (t) + \H(t) is a solution of the differential equation 

C"{t) = -47T 2 H(t)-7TH'{t). 

1 

Corollary: A§(0) = /i 3 (0) + 2^ 2 / H 2 (t)dt. 

o 

In particular, for all parameters E ^ near zero we Ziaue £/ie inequality 

fi 3 (E) < X 2 (E). 

Moreover, the first positive eigenvalue Hi(gE) of the Laplace operator is always 
smaller then the corresponding eigenvalue X\{gE) of the Dirac operator for any metric 
gE near E « 0, i.e., 

Hi(9e) < A?(5b)- 

Remark: The explicit formulas in Theorem 1 and Theorem 2 can be generalized 
to the case of an arbitrary conformal deformation. Fix a Riemannian metric g Q on a 
surface M 2 and consider the deformation 

g E = (1 + EH + E 2 G + 0(£ 3 ))<7o 

of the metric. Moreover, suppose that fii(E) is the deformation of the eigenvalue of 
the Laplace operator and is the corresponding family of eigenfunctions. Then the 
following formulas hold: 

J Hf 2 dM 2 

a. ) £i(0) = -Mi(0) 

M 2 

J (Gf 2 + Cf )dM 2 

b. ) Ai(0) = -2 M i(0) , 

/ fodM 2 

M 2 

where the function C is the solution of the differential equation 

A o C = /x 1 (0)iJ/ o + /x 1 (0)C. 

The corresponding expression for the variation of the eigenvalue Xi(E) of the Dirac 
operator can also be computed: 
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I H-\i) \ 2 dM 2 

c.) Ai(0) = -Ai(0) 

2 J \^ \ 2 dM 2 

M 2 

and a similar formula holds for the second variation. 

Once again, a similar, though even more intricate computation yields the fourth vari- 
ation of X 2 (E) under the assumption that all previous variations of X 2 (E) vanish. 
This is needed for the discussion of the example in Section 5. 

Theorem 3: Consider a deformation 

9E = {1 + EH(t))g 

of the flat metric g Q on the torus T 2 and suppose that the following conditions are 
satisfied: 

a. ) H(t) = H(l-t); 

b. ) A|(0) = A|(0) =A|(0)=0. 

Then the fourth derivative [A§(0)] (IV) of the spectral function X^(E) at E = is given 
by the formula 

l l l 

[A|(0)] (IV) = Q J H 3 (t)H"(t)dt + y J H 2 (t)(H'(t)) 2 dt - J (l6Tr 2 H(t) + H"(tfjC 3 (t)dt 



1 

+ / Q(H'(t)) 2 + 2H(t)H"(t) ) jC 2 (t)dt 
o 

i 

-J (l5(H'(t)) 2 +6H"(t)H(t))H(t)C 1 (t)dt, 
o 

where the functions Ci(t),C2{t),C 3 (t) are periodic solutions of the equations 
C'((t) = -4ir 2 H(t)-^H"(t)-7rH'{t) 

C 2 \t) = \H(t)H"{t) + b -{H'(t)) 2 + 2irH\t)H{t) - (Vtf (t) + l -H"{t) + 27rff'(*)) d(t) 
Cm = - ^-H"(t)H\t) + l ^H(t)(H'(t)) 2 + ^H'(t)H 2 (t)) 

+ (^H(t)H"(t) + l(H'(t)) 2 + QH{t)H'{t)) C 1 (t) 

- (3*H'(t) + hl"(t) + 47r 2 F(t)) C 2 (t). 
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Remark: In the special case of H"(t) = — 16ir 2 H(t) the derivative [A 2 (0)] (IV) does 
not depend on Cs(t) and the formulas become much simpler. Such a metric will be 
the object of Section 5. 



4 Examples 

4.1 The variation g E = (1 + E cos(2irt))g 

The volume vol (T 2 ,gE) = 1 of this variation of the flat metric g Q is constant and all 
first derivatives at E = vanish since 

l l 

J cos(27rt)cos 2 (2vrt) = J cos(2vrt) sin 2 (2vrt) = 0. 


A computation of the second derivatives yields the following numerical values: 

2 10 
/ii(0) = --vr 2 , il 2 (0) = -7r 2 , /i 3 = -47r 2 

A?(0)=vr 2 , A 2 (0) = A 2 (0) = -3n 2 . 
In particular, we obtain 

for all parameters i? / near zero. The eigenspinor corresponding to the minimal 
positive eigenvalue of the Dirac operator does not vanish anywhere (Figure 1). 




(Figure 1) 
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4.2 The Mathieu deformation g E = (1 + E cos(4irt))g of the fiat 
metric 

This deformation of the flat metric again preserves the volume, and the Laplace 
equation essentially reduces to the classical Mathieu equation u"(x) + (a + 
16g cos(2x))u(x) = 0. In this case the first variation is trivial only for the Dirac 
equation. Indeed, we have 

/i 1 (0) = 2vr 2 , / l2 (0) = -2vr 2 , /i 3 (0) = 

Af (0) = A 2 (0) = A|(0) = 0. 
Even for the Mathieu deformation we conclude that 

Hi{9e) < ^i(9e) 

for E 7^ near zero. A computation of the second derivatives yields the following 
result (Figure 2): 

/i 3 (0) = -vr 2 , A 2 (0)=vr 2 , A 2 (0) = A|(0) = 0. 
For a detailed discussion of this metric, we refer to the next section. 



U2(E) m(E) 




(Figure 2) 
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4.3 The variation g E = (1 + E cos (2tt Nt)) g , N > 3 

Since 

l l 
J cos(2ir Nt) cos 2 (2irt)dt = J cos(2vriVt) sin 2 (2vrt)dt = 



for N > 3, the first variations of our spectral function vanish. We compute the 
second variation using the algorithm in Theorem 2: 

4vr 2 4vr 2 
j il (p) = j i2 (p) = - w — , /i 3 (0) = - — 

A 2 (0) = vr 2 , A 2 (£) = A|(0) = (l - A) vr 2 . 
In particular, we obtain again 

>?i(9e) > Mi(£te) 
for all parameters E ^ near zero (Figure 3). 





*1(E) 
x i 39.49 

X \ 
X \ 

X \ 

X V 
X^ X 

\ ^39.485 

x \ 

\ 

39^48 

NX 


X2(E)=^3(E) 

i i 
i i 
i i 
i i 

1/ 

- / / 




//S9.475 




0.1 


/-0/05 


0\05\ 0.1 



(i 3 (E) Hi(E)=n 2 (E) 



(Figure 3) 
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5 The Mathieu deformation of the flat metric 



In the previous examples, the deformation 

g E = (1 + E cos(4wt))g 

of the fiat metric g Q plays an exceptional role, because the derivatives /ii(0), /i2(0) ^ 
are non-zero. Therefore, we study the behaviour of the first positive eigenvalue for 
the Laplace and Dirac operator in more detail. First of all, the lower bound 

4tt 2 

' < Vi(E), Ai(-E) 



/i 4 



max 



yields the estimate 



4tt 2 

" <^(E\\l{E) 



l + \E\ 



for all parameters — 1 < E < 1. In case of the function f(t) = sin(27rt) the upper 
bound B£(g, f) of Section 2 leads to the estimate 



2 + \E\' 

i.e., for all parameters — 1 < E < 1 the inequality 

4vr 2 ,^ 8vr 2 

< Mi(^) < 



l + \E\ - ^ y ' ~ 2 + |£| 
holds. On the other hand, for the Dirac operator the function 

fit) = (1 + £cos(47rf))2 sin(27rt) 
gives an upper bound Bf)(gE,f) for its first eigenvalue with the property 

\im B u D (gE,f)=5n 2 . 
h/ — ► — 1 



We will thus investigate the limits lim ui (E) as well as lim \%(E). The eigenvalue 

E->-l E->-l 

fJ-i(E) is related with a periodic solution of the Sturm-Liouville equation 
A"(t) = —/j,i(E)(l + £'cos(47rt)^y4(t) + 4ir 2 k 2 A(t), 

where k = 0, ±1 (see Proposition 1). Let us introduce the function B(x) := A (^%^ 
where < x < 2ir. Then the Sturm-Liouville equation is equivalent to the classical 
Mathieu equation 

B"(x) + (a + 16qcos(2x))B(x) = 0, 
where the parameters a and q are given by 
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a " 4vr2 fc ' 9 "16(4^) ' fc " U ' ±L 

For E — > — 1 the parameters of the Mathieu equation are related by 

o = -16<7-/c 2 , fc = 0,±l. 

Using the estimates for ni(E) we obtain 

2vr 2 < lim u^E) < V, 

hi — > — 1 o 

i.e., — < n < _ — i n case E = —1. 

24 _ H ~ 32 

A numerical computation shows that, under these restrictions, the Mathieu equation 
has a unique periodic solution for k = and q ~ 0,04113. This solution B(x) is 
the first Mathieu function se±(x, q), which is the deformation of the function sin(x). 
Consequently, we have 

lim m(E) = -16 • q ■ 4vr 2 « 2, 6323vr 2 . 

— > — 1 

The limits of the spectral functions /^(-E) and ^{E) can be computed in a similar 
way: 

lim (i2(E) ~ 1, 79 ■ (4tt 2 ) , lim /x 3 (£) « 0, 9 • (4^ 2 ). 

E — > — 1 E — > — 1 

These limits correspond to the Mathieu functions ce\(x,q) and ce Q (x,q) for the pa- 
rameters q Ri —0, 112 in case of /-12(E) and q w —0, 056296 in case of [iz(E). 




(Figure 4) 
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Approximation of the periodic solution for [ii(E),E — ► —1 : 

NDSolveKy" [x] + 32 (0 . 04113) (Sin [x] ) "1 y[x] == 0, 
y[x] == , y> [0] == 1} , y , {x , , 10 Pi}] 

Plot [Evaluate [y[x]/.y„ , {x , , 10 Pi}] 




(Figure 5) 



Approximation of the periodic solution for /j,2(E),E — ► —1 : 

NDSolveKy" [x] + 32(0. 1112) (Sin [x] ) "2 y[x] == 0, 
y[x] == 1 , y> [0] == 0} , y , {x , , 10 Pi}] 

Plot [Evaluate [y[x]/. % , {x , , 10 Pi}] 



5 
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(Figure 6) 



Approximation of the periodic solution for ^{E),E 



-1 



NDSolveKy" [x] + 32 (0 . 056296) (Sin [x] ) "2 - 1) y[x] == 0, 
y[0] == 1 , y> [0] == 0} , y , {x , , 10 Pi}] 

Plot [Evaluate [y[x]/.% , {x , , 10 Pi}] 

A A A A A A A A 

1.5 



1 . 4 
1 . 3 
1.2 
1 . 1 



10 15 20 

(Figure 7) 
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The eigenvalues \ 2 a {E) of the Dirac operator are related with the periodic solutions 
of the Sturm-Liouville equation 



-A"(t) = \ 2 h\t) + 



h(t)h"(t)-2(h'(t)) 2 A 2l2 , A h'(t) 



h 2 (t) 



4tt 2 1 2 +4ttI-^\ A(t). 
n(t) I 



For the Mathieu deformation we have 



h{t)h"{t) - 2{ h'(t)) 2 



-4it 2 E 



E + cos (47rt) + |£sin 2 (47rf) 
(l + ^cos(47rt)) 2 ' 



First we discuss the case that I = 0. Then the first positive eigenvalue of the Dirac 
equation is given by 



A 2 = 



4tt 2 



2 ■ 



h 2 (t)dt 



\0 



In case of the Mathieu deformation we obtain 

l l 



2^2 



^lin^y h 2 (t)dt = J yjl - cos(47rt)dt = 
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and, finally, 

£ lim i A 2 (i?) = V^(4,92)vr 2 




lower bound 

-0.~6 -0.4' 

(Figure 8) 



4 7I 2 



We now investigate the case 1 = 1. Let us consider the Hamiltonian operator He 
given by the Sturm-Liouville equation for A 2 = 0: 



H E 



dt 2 



+ 47r 2_ 4n h'(t) h{t)h"{t)-2{h'{t)f 



dt 2 



+ PE(t), 



h(t) h 2 {t) 

where the potential PE(t) is given by the formula 

2 4cos(4vrf) + 4sin(4vrt) + £sin 2 (4vrf) + 2E(2 + sm(8vrf)) 



p E (t) = Air 2 + Eir' 



(1 + £cos(47rf)) 2 



For all parameters — 1 < E < the Hamiltonian operator He is strictly positive (see 
Proposition 3). Consequently, the eigenvalue A3 (£7) is the first number such that 

inf spec {H E - A 2 (l + £cos(47rf))) = 0, 
and the corresponding solution of the Sturm-Liouville equation 

A E (t) = (PE(t) - Xt(E)(l+Ecos(ATTt)))A E {t) 
is unique and everywhere positive. In particular, the solution satisfies the condition 

A E {t + \) = A E {t). 
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Since ^4_e(t) is a positive periodic solution of the Sturm-Liouville equation, we obtain 
the condition 

1 

J (PE(t) - A§(£7)(1 + E cos(4irt)))dt > 
o 

and thus an upper bound for X^(E): 

l 

Xl(E) < J PE (t)dt. 
o 
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40.5 








40.25 




-0 . 5 


-0.4 -0\3 -0.2 


-0 . 1 








39.75 





(Figure 9) 



We notice that this upper bound for X^(E) grows and reflects, indeed, the real 
behaviour of X^(E) near E = 0. To see this, we use Theorem 3 to compute the 
fourth variation of this spectral function (the third variation vanishes since X\(E) 
has to be a symmetric function in E). One obtains the following result: 

[A|(0)r> = ^vr 2 > 0. 

On the other hand, using well-known approximation techniques for Sturm-Liouville 
equations with periodic coefficients (see [9]) we can approximate X^(E) for a fixed 
parameter E. Indeed, one replaces the potential in the Sturm-Liouville equation by 
the first terms of its Fourier series. This reduces the computation of the approxima- 
tive eigenvalue to a finite-dimensional eigenvalue problem. For example, in case of 
E = —0.3 the mentioned methods yields the result 

A|(-0.3) » 39.6733. 

Let us study the behaviour of the spectral function X^(E) for E — > —1. More 
generally, denote by X 2 (E,l) the first eigenvalue of the Dirac operator such that the 
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corresponding eigenspace contains an S^-representation of weight I. In particular, 
we have X 2 (E) = X 2 (E, 1) = X 2 (E, — 1). We apply the Corollary of Proposition 3 to 
the function h^it) = \J\ + E cos(47rt) and conclude that 

/ ^!g2^£H - A»(B,I) / Vl + Bccf^f.)* > 
J Jl + £cos(47ri) 7 v 

v v y 

holds for any periodic function ip(t). Fix a test function ip(t) and consider the limit 
E — ► — 1. Then we obtain the inequality 

1 (27r^(t) - y/(t)) 2 
|sm(2vrf)| 



lim A 2 (£,Z) < i ^- r 



J \sm(2irt)\(p 2 (t)dt 
o 

We apply this estimate to the function 

. . cos(27rt) + I sin(27rf) 

= — 27pTI)^ — • 

Then 2irl<pi(t) — f'i(t) = sin(27rt) and we obtain the following 
Proposition 4: 

Remark: At E = we have A 2 (0,Z) = 47r 2 / 2 . On the other hand, for / > 3 the 
inequality 

1 + 2P 

holds, i.e., 



limA^E,*) < A 2 (0,/) Z>3. 

The latter inequality means that the eigenvalue X 2 (E,l) decreases for E — > — 1 
(i > 3). 

The behaviour of X 2 (E) = X 2 (E,l) for I = 1 is completely different. This spectral 
function increases for i£ — ► — 1. Using the formula 

1 (27TM*) " ^)) 2 



X 2 3 (E) = inf ^ 



^1 + £cOs(4vrt) 



y y / l + J Bcos(47rt)v3 2 
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we can approximate the positive minimizing Mathieu spinor MS (E, t) of topological 
index I = 1 by expanding it in its Fourier series. We thus obtain for example: 

E = -0.9: A§(-0.9) « 40.1464 

MS(-0.9,t)=(Sqrt[Sqrt[l+ (-0.9) Cos [4 Pi t] ] ] ) Sqrt [1+ (0.44)Sin[4 Pi t] 
+ (0.15)Cos[4 Pi t] + (0.09)Sin[8 Pi t] + (0.17)Cos[8 Pi t] 
+ (0.028)Sin[16 Pi t] + (0 . 051)Cos [16 Pi t] + (0 . 051) Sin [12 Pi t] 
+ (0.085)Cos[12 Pi t] + (0.016)Sin[20 Pi t] + (0 . 026) Cos [20 Pi t] 
+ (0.01)Sin[24 Pi t] + (0 . 014)Cos [24 Pi t] + (0 . 005) Sin [28 Pi t] 
+ (0.007)Cos[28 Pi t] + (0 . 0033) Sin [32 Pi t] + (0 . 0044) Cos [32 Pi t]] 

E = -0.95: A§(-0.9) « 44.6024 

MS(-0.95,t)=(Sqrt[Sqrt[l+ (-0.95)Cos[4 Pi t]]])Sqrt[l+ (0.585)Sin[4 Pi t] 

+ (0.049)Cos[4 Pi t] + ( 0.1)Sin[8 Pi t] + (0.08)Cos[8 Pi t] 

+ (0.063)Sin[12 Pi t] + (0 . 063)Cos [12 Pi t] + (0.041)Sin[16 Pi t] 

+ (0.04)Cos[16 Pi t] + (0.026)Sin[20 Pi t] + (0 . 026)Cos [20 Pi t] 

+ (0.017)Sin[24 Pi t] + (0 . 018)Cos [24 Pi t] + (0 . 012) Sin [28 Pi t] 

+ (0.011)Cos[28 Pi t] + (0.008)Sin[32 Pi t] + (0 . 007)Cos [32 Pi t]] 



Finally, we can compute the limit lim X^(E) replacing again the potential in the 

E — > — 1 

Sturm-Liouville equation by the first terms of its Fourier series. For E = — 1 this 
amounts to studying the differential equation 

sin 2 (27rt)^"(t) = |t^ 2 (9 - 3cos(4vrt) - 4sin(4vrt)) - 2A§ sin 4 (2vrf)J A(t) 

and the finite-dimensional approximation yields the result 

lim Xl(E) « 47.2437. 

E — > — 1 

Remark: The second variation formulas prove that, in case of the family 
qe = (l + -E , cos(27rt))<7 (N = 1), the minimal positive eigenvalues of the Laplace and 
Dirac operator decrease (see Example 4.1) and are smaller than 47r 2 . The numerical 
evaluation of fJ^(E) and X^(E) yields the following table: 



E 





-0.1 


-0.3 


-0.5 


-0.7 


-0.9 


-0.95 


-0.99 


-1 




4vr 2 


39.284 


37.897 


35.741 


33.378 


31.09 


30.5 


30.1 


30.013 


Ai 


4tt 2 


39.333 


38.353 


36.714 


34.983 


33.331 


33.2830 


36.04 


« 36.2 
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6 Final remarks 



As shown previously, any local deformation gE of the flat metric realizes the inequal- 
ity 

Hi(9e) < Ai(flfc) 

between the first eigenvalues of the Laplace and Dirac operator up to second order. 
We are not able to give an example of a Riemannian metric g on T 2 such that 
\i(g) < fJ-i(g) holds. Moreover, denote again by Af(p;Z) the first positive eigenvalue 
of the Dirac operator such that the eigenspace contains an S^-representation of weight 
/ € Z. The corresponding eigenvalue of the Laplace operator we shall denote by 
fj,i(g;l). It is a matter of fact that in all families of Riemannian metrics we have 
discussed these two eigenvalues are very close. Let us consider, for example, the 
metric gE by the function 

h E (t) = e? (sin(27rt)-2cos(27rt))_ 

For the parameter E = 1 we obtain the following numerical values using the 
approximation method described before in the space spanned by the functions 
1, sin(2-7rnt), cos(2-7rnt) (1 < n < 5): 

\\{g]l) « 6.11056 , m{g;l) » 5.19025. 

However, even in this case we already have the inequality Hi(gE', 1) < Af (gE', 1) and 
the following figure shows the graph of the two spectral functions for < E < 1 (for 
the first and the second positive eigenvalue): 




(Figure 10) 
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